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ABSTRACT 

Continuing our study of the role and evolution of accretion disc winds in the micro- 
quasar GRS 1915+105, we present high-resolution spectral variability analysis of the 
/3 and 7 states with the Chandra High Energy Transmission Grating Spectrometer. 
By tracking changes in the absorption lines from the accretion disc wind, we find new 
evidence that radiation links the inner and outer accretion discs on a range of time- 
scales. As the central X-ray flux rises during the high-luminosity 7 state, we observe 
the progressive over-ionization of the wind. In the (3 state, we argue that changes in 
the inner disc leading to the ejection of a transient 'baby jet' also quench the highly- 
ionized wind from the outer disc. Our analysis reveals how the state, structure, and 
X-ray luminosity of the inner accretion disc all conspire to drive the formation and 
variability of highly-ionized accretion disc winds. 

Key words: accretion, accretion discs - stars: winds, outflows - black hole physics - 
instabilities - X-rays: binaries - X-rays: individual (GRS 1915+105) 
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1 INTRODUCTION 

In the last two decades, ASCA, Chandra, XMM-Newton, 
and Suzaku have discovered a multitude of highly ionized 
absorbers in the X-ray spectra of black hole and neutron 



for b lue shifts (e. 



20091 ). 



In some cases, without clear evidence 


Ebisawa 19971: 


Kotani et al.l 19971. 


20021; iMiller et al. 


2004 Xiane: et al.l 



these absorbers have been attributed to hot or 
extended disc atmospheres. Many observations, however, 
provide unambig uous evidence for outflo w s from the accre- 
tion d isc (e.g. iBrandt fc Schulj l2000l: ISchulz fc Brandt] 
20021: lUedaetalJ 12004 IMiller et al l l2006allbl. 12008 



Ueda. Yamaoka. fc Remillardl 120091; INeilsen fc Led 12009 
Blum et all |2010|; iRevnolds fc Millerl l20ld ; IMiller et all 
201ll : INeilsen, Remillard. fc Led l201ll . hereafter Paper II). 



These results suggest that accretion disc winds may be as 
common as relativistic jets in accreting systems. 

Of particular interest, then, is the variability of these 
disc winds and their role in the evolution of the accretion 
flow. Ionized absorbers in X-ray binaries are now known 
to vary on time-scales from 5 seconds (GRS 191 5+105, 
iPaper ij) to 300 s (H 1743-322. IMiller et all l2006tj ) to ks 
or more (Cir X - l, ISchulz fc Brandtl |2002t GRS 1915+105, 
iLee et all 120021 : lUeda et al.l |2010T ). Furthermore, there is 
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mounting evidence that this variability cannot be due 
solely to changes in the ionizing luminosity (ILee et al.ll200^ ; 
ISchulz fc Brandtl |2002| ; iBlum et all l20ld ; IPaper it ), which 
implies that high-spectral-resolution X-ray observations of 
accreting X-ray binaries may be able to probe either the 
physical processes that link outflows to the behavior of the 
inner accretion flow, or inhomogeneities in the structure and 
density of accretion disc winds. 

In a study of 10 years of Chandra Hig h-Energy Trans- 
missi on Grating Spectrometer (HETGS; ICanizares et alj 
120051 ) observations of GRS 1915+105, we demonstrated 
that the strength of absorption lines in the accretion disc 
wind is anticorrelated with the fractional hard X-ra y flux 
(and the stre ngth of the jet; INeilsen fc Led 120091 . here- 
after IPaper J ; for a preliminary discussion of the absorp- 
tion lines, see IMiller et al.ll2008f ). Based on this anticorrela- 
tion and estimates of the mass loss rates in the wind and 
the jet, we argued that accretion disc winds may be able 
to suppress or quench jets by directly draining their mat- 
ter supply on long time-scales. The fractional hard X-ray 
flux is therefore a useful indicator of both the accretion 
state and the physics of outflows around black holes. Yet 
much more can be said for GRS 1915+105, whose X-ray 
lightcurve exhibits a t least 14 differe nt patterns of variability 
dBelloni et all |2000|. her eafter BOO; iKlein-Wolt et all 120021 ; 
lHannikainen et al.l 120051 ) . Many of these phenomenological 
patterns, which are highly structured and high- amplitude, 
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Figure 1. The Chandra (black) and RXTE (red) lightcurves of the (top, middle) and 7 (bottom) states. The RXTE lightcurves have 
been renormalized to match the mean HETGS count rate, with excellent agreement between the two missions. In the top and middle 
panels, we use horizontal lines to indicate the locations of the X-ray dips analysed in Section l3.ll 



are interpreted as limit cycles of accretion and ejection in 
an unstable disc jBelloni et al1ll997l: iMirabel et all 1 19981 : 
iTagger et ail 12004 iFender fc Bellonil I2004T ). In order to un- 
derstand the physics linking disc winds, jets, and this fast X- 
ray variability, we presented the first detailed phase-resolved 
spectral analysis of a variability class in GRS 1915+105, us- 
ing a joint RX TE/Chandra HETGS observation of the p 
state ( Paper IIP . We showed for the first time that changes 
in the broadband X-ray spectrum produce changes in the 
structure and mass loss rate of the accretion disc wind on 
time-scales as short as 5 seconds. Furthermore, our anal- 
ysis indicated that the wind in the p state may be suffi- 
ciently mas sive to excite long-term oscillations in t he mass 
of the disc (|Shields et all 1 19861 ; iLuketic et a"i1l20101 l. By es- 
tablishing clear ca usal links between the X-ray oscillations 
and the disc wind, iPaper ij provides a physical mechanism 
that could be responsible for both the suppression of jets 
and (possibly) state transitions in black hole binaries. 

In this paper, we continue our follow-up wind variabil- 
ity studies with spectral timing analysis of three additional 
archival Chandra HETGS observations of GRS 1915+105. 
Two of these observations (ObsIDs 6579 and 6580) fall into 
the (3 class, a wild ~ 30-min cycle that is associated with 



Table 1. Chandra HETGS observations of GRS 1915+105 



(iPoolev & Fendei 1997 


Fender et al.l 19971; Mirabel et al.l 


1998; Eikenberrv et alj 


Il99gj; IMirabel & Rodriguez! |l999l; 


Belloni et al. 2000; Klein- Wolt et al. 2002: Fender & Bellonil 



Obs. 


X-ray 




Elapsed 


Tcxp 


ID 


State 


Date 


Time (ks) 


(ks) 


6579 


13 


2005 Dec. 1 


20.65 


12.30 


6580 




2005 Dec. 1 


21.62 


12.14 


6581 


7 


2005 Dec. 3 


20.65 


9.73 



2004). Each P cycle begins with a 10-15 minute dip 
in the lightcurve, which is spectrally hard and exhibits 
a QPO whose frequency tracks the X-ray flux (e.g. 
Markwardt et aJj|l999; iMikles et all 120061 ; iRodriguez et al.l 



2008bh . The dip ends with a bright spike, followed by 10- 
15 minutes of soft, bright X-ray flares. These flares are 
coincident with strong optically-thin radio/infrared oscil- 
lations (baby jets), but there is an ongoing debate as to 
whether these discrete ejections are produc ed during the 
dip i ts elf or at the time of the X-ray spike jMirabel et alj 



1998; 



Eikenberrv et ail Il998l: iMig lari fc Bellonil I2OO3F 
Rothstein ct al . 2005]; IRodriguez et alJl2008al ) 

During the third observation (ObsID 6581), GRS 
1915+105 exhibited irregular variability with very high X- 
ray flux (RXTE count rate > 6000 counts s _1 PCU" 1 ) and 
extremely strong red noise in the power density spectrum 
( P Y ~ u~ 2 ). Although we identified this as the 7 state 
in IPaper J and continue to refer to it thus here, this par- 
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ticular observation is rather unusual since it also displays 
significant similarities to typical S and fi states. Relatively 
little has been said about the /j,, 5, and 7 states in the 
liter ature, as none have bee n associated with strong radio 
jets dKlein-Wolt et al]l2002l'). However , in a study of the fj, 



Soleri. Belloni. fc Casellal (|2008h reported the discov- 



state, 

ery of a transient low-frequency QPO, which they argue may 
be linked to state transitions or transient jets. This leaves 
open the possibility that the /x/<5/7 states exhibit physically- 
interesting outflow behavior, despite their weak radio flux. 
In this respect, these states could b e similar to the p state , 
which has very low radio brightness (|Klein-Wolt et al]l2002l 'l 
but exhibits ext reme disc wind variability and may produce 
short-lived jets (|Paper 111) . 

Here we address the variability of the known accretion 
disc wind in the Chandra HETGS observations of the /3 and 
7 states, whose lightcurves are shown in Figure [T] By com- 
paring changes in the ionized absorption to changes in the 
X-ray flux, we show how the wind is linked to the oscilla- 
tions in the inner disc. The paper is organized as follows: in 
Section [5] we describe our observations and data reduction. 
In Section [3] we perform spectral analysis of the wind in the 
P state with the Chandra HETGS and RXTE PCA; ing]we 
use the HETGS to study the wind in the 7 state. We discuss 
the implications of our results in Section [5] 

2 OSERVATIONS AND DATA REDUCTION 

GRS 1915+105 was observed with the Chandra HETGS on 
2005 December 1 (01:45:09 UT and 18:41:30 UT) and Dec 

3 (16:26:06 UT) for 20.65 ks, 21.62 ks, and 20.65 ks, respec- 
tively. The data were taken in Continuous Clocking Mode, 
which has a time resolution of 2.85 ms, in order to mitigate 
photon pileup. However, due to the high X-ray flux the ob- 
servations suffered from severe telemetry saturation, so that 
the good exposure times for these observations were only 
12.30 ks, 12.14 ks, and 9.73 ks (details in Table [TJ. 

We reduce and barycentre-correct the Chandra data us- 
ing standard tools from the CIAO analysis suite, version 4.3. 
We use the order-sorting routine to remove the ACIS S4 
readout streak, since the destreak tool can introduce spectral 
artefacts for bright continuum sources like GRS 1915+105. 
After reprocessing and filtering, we extract High-Energy 
Grating (HEG) spectra and create grating responses. For 
the 7 observation (ObsID 6581), we find that the detected 
position of the th -order image results in an offset of —44 eV 
between the HEG —1 and HEG +1 orders (measured from 
the centroid of the Fexxvi line at ~ 7 keV). Shifting the po- 
sition of the th -order image to the nominal source position 
results in an offset of 36 eV, so we take the average of the de- 
tected position and the nominal position as the true location 
of the th -order image on the CCD. The resulting position 
error from this average is less than 0.1 pix, so the associ- 
ated uncertainties in the wavelength calibration are smaller 
than our uncertainty in the line centroids. We also extract 
10-second lightcurves with dmextract. 

Due to incomplete calibration of Charge- Transfer In- 
efficiency (CTI) in CC mode, there is some wavelength- 
dependent disagreement in the continuum flux between spec- 
tral orders of the HEG (and the MEG, which we do not con- 
sider here because of its lower spectral resolution). For this 
reason, it is not currently possible to fit a physical contin- 
uum model to the HETGS data. Instead, we fit the individ- 



ual spectra with polynomials to model the local continuum 
and use Gaussians for line features found in the combined 
residuals. 

During the Chandra observations, RXTE made pointed 
observations of GRS 1915+105, lasting 23.1 ks, 13.7 ks, and 
21.7 ks (elapsed), with exposure times of 13.2 ks, 7.2 ks, and 
13.6 ks, respectively. We select all available data subject to 
the following constraints: (1) the Earth-limb elevation angle 
is above 3°; (2) the spacecraft is outside the South Atlantic 
Anomaly; (3) the offset angle from GRS 1915+105 is less 
than 0.02°. In this paper, we analyse the data from the Pro- 
portional Counter Array (PCA), which covers the 2-60 keV 
band, in order to verify and supplement the variability de- 
tected by Chandra. For timing analysis, we make use of the 
data from the binned mode B_8ms_16A_0_35_H_4P, which 
covers the 2.0-14.8 keV band at 7.8 ms time resolution, and 
the event mode E_16us_16B_36_ls, which covers the 14.8-60 
keV band at 15.3 fis time resolution. We extract and com- 
bine 1-second background subtracted lightcurves from each 
of these modes. For comparison to the Chandra data, we 
barycentre the RXTE lightcurves. For spectral analysis, we 
extract Standard-2 129-channel spectra from relevant time 
intervals (see Section and Figure [TJ. We restrict our anal- 
ysis to the 3-45 keV top layer spectrum of PCU-2, which 
is best calibrated. We assume 0.6% systematic errors. Al- 
though these states exhibit some strong variability on time- 
scales shorter than the 16-second bins of the Standard-2 
data, our goal is to use the spectra to understand the ion- 
ization of the wind, not to characterize physical variability 
in the accretion disc and corona. 

3 THE STATE 

In this section, we explore the variability and dynamical evo- 
lution of the accretion disc wind in the P state. Here, links 
between the extreme X-ray variability and the disc wind are 
particularly interesting given the association of the P state 
with transient jets and disc-jet interactions. Previously, we 
analysed the average spectrum of these two (ObsIDs 6579 
and 6580; Table [T]) joint RXTE/Chandra observations of 
GRS 1915+105 in the P state (Paper I). We found the av- 
erage RXTE continuum during each of these observations 
to be relatively soft, with ~ 79% of the 3-18 keV X-ray lu- 
minosity emitted below 8.6 keV (Lx ~ 6 x 10 38 ergs s _1 ). 
In both time-averaged high-resolution X-ray spectra from 
the HETGS, we detected an Fexxvi Lya absorption line 
from the accretion disc wind. This feature had equivalent 
widths of W = 13.3±l g eV and 19.3±|-§ eV in observations 
6579 and 6580, respectively; the corresponding blueshifts 
were 910±^° km s _1 and llOOtggo km s_1 - In what fo1 " 
lows, we use the known properties of the P state to extract 
high-resolution X-ray spectra corresponding to physically- 
interesting intervals of this unusual oscillation. We com- 
pensate for severe telemetry saturation (leading to reduced 
counts, Section [2] Table [T]) by fitting both observations 6579 
and 658 jointly. All spectral fit t ing is done f rom 5-8 keV 
in ISIS (|Houck fc Denicolal l200d : iHouckl 12002] ). We assume 
a distance and inc lination of D = 11.2 kpc an d i = 66° 
dFender et alJll999h : we fix N H = 5 x 10 22 cm" 2 (|Lee et al.l 
120021 and references therein). 

In general, choosing time intervals for spectral timing 
analysis is difficult, but the P lightcurve is easily divided 
into two sections with different physical properties: (1) the 
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Table 2. X-ray absorption lines in GRS 1915+105 



Eq E obs At; shift <7 V Wo 



Interval 


F 5 _g 


Line 


(keV) 


(keV) 


(km s 1 ) 


(km s 1 ) 


Flux 


(eV) 


Dip 


6.0 ± 0.5 


Fe xxvi 


6 


.961 




-1340±gg 


< 1310 


1 n+0.7 
L9 -0.6 


1 1 Q + 4.1 






Fexxv 


6 


.700 


6.730^ 


— 1340 a 


= Fexxvi 


< 0.8 


< 4.8 


Flare 


7 fi _i_ n /i 
/ .0 it U.4 


Fe xxvi 


6 


.961 


fi QQ9+0-006 

o.y»2_ 005 


onn+230 


fifin + 230 

6bU_ 240 


D.4 it U.o 


ofi n+3.8 
2b.U_ 3 g 






Fe xxv 


6 


.700 


c 71 o+O.OlO 
O.' 18-0.008 


rron + 360 

-780I 450 


= Fexxvi 


2.7 it 0.7 


i 1 zi+2.8 


Pi 


A fi + A 


Fe xxvi 


6 


.961 


D-y»O_0.010 


1060+ 420 




^ 9 -t- n i 

iJ.Z It U. 1 


AO 9 + 8-1 






Fexxv 


6 


.700 


6.71 it 0.01 


- 560 -590 


= Fexxvi 


3.5+'^ 


24.6+°° 


Pi 


0. 1 III u.o 


UVi W \ / T 

re xxvi 


6 


.961 


fi Q84+ - 012 


-yt>u_ 520 


S 700 


Qoj.ni 

O.O It U. 1 


99 n + 5.9 






Fexxv 


6 


.700 


„ r71 o +0.013 


or.n+390 
- 800 -560 


= Fexxvi 


2.1 it 0.9 


t 1 o+4.8 


Q 

P3 


lo.y ± 0.3 


Fe xxvi 


6 


.961 


fi oan+0-008 
6.980_ 005 


-'° U -360 


< 910 


8.6_ 2 .0 


r,r> o+4.9 

22 - 3 -5.1 






Fexxv 


6 


.700 


6.718° 


-780 a 


= Fexxvi 


o n+1.7 
° -1.8 


7 c+3.9 
'■°-4.1 


71 


18.8 it 1.5 


Fe xxvi 


6 


.961 


6.985 it 0.008 


-I020±lf 


■7nn+410 
700 - 430 


11.6 it 2.8 


21.3 it 5.1 






Fexxv 


6 


.700 


6.723 a 


= Fexxvi 


= Fexxvi 


q k+2.2 
0.0_ 2 3 


5 q+ 3 - 7 


72 


22.5 it 2.0 


Fe xxvi 


6 


.961 




-9201^0 


980t|« 


«-l±3.0 








Fexxv 


6 


.700 


6.721 a 


= Fexxvi 


= Fexxvi 


< 5.7 


< 7.9 


73 


24.8 it 2.2 


Fe xxvi 


6 


.961 


6.990 it 0.007 


-1240 it 300 






26.4+5-1 






Fexxv 


6 


.700 


6.728° 


= Fexxvi 


= Fexxvi 


< 4.5 


< 5.7 


74 


26.4 it 2.3 


Fe xxvi 


6 


.961 


6.998 it 0.007 


— 1550^280 


7gQ +490 
'» u _680 


17.6+ 3 3 ;? 


oq i+5.1 
zc> - ± -4.8 






Fexxv 


6 


.700 


6.735 a 


= Fexxvi 


= Fexxvi 


< 6.6 


< 7.6 


75 


28.3 it 2.2 


Fe xxvi 


6 


.961 


6 gg6+ 008 


-1490+ 440 


< 1330 


I2.9l§;5 


i fi n+ 3 ' 5 






Fexxv 


6 


.700 


6.734 a 


= Fexxvi 


= Fexxvi 


< 6.5 


< 7.4 



Errors quoted are 90% confidence ranges for a single parameter unless otherwise stated. State: X-ray variability type in 
the classification of iBOOl . Interval: line parameters are reported for different flux levels and time intervals (see text for 
more details). -F5-8 : 5-8 keV (absorbed) continuum flux in units of 10 — 9 photons s — 1 cm -2 , with errors given by the 
standard deviation between the HEG ±1 orders; Eo : rest energy; -E bs ! measured energy; Ai> s hift : measured Doppler 
velocity; <t v : line width; Flux: measured absorbed line flux in units of 10 — 3 photons s _1 cm -2 ; Wo '■ line equivalent 
width. 

a Since the Fexxv line is not very significant during this interval, we tie its velocity to that of Fexxvi. 



long, hard, discrete jet-producing dip and (2) the subsequent 
period of bright flaring. Examination of the /3-state X-ray 
lightcurves in the literature and in our data indicates that 
the hard dip is bracketed on the left by the last bright flare 
of the preceding flaring interval, and on the right by peak or 
falling edge of the X-ray spike. Using the Chandra lightcurve 
and the PCA lightcurve where available, we hand-select the 
dip and flare intervals accordingly. The resulting GTIs are 
shown underneath the X-ray lightcurves in the top and mid- 
dle panels of Figure [T] 

3.1 Chandra HETGS 

Applying these GTIs to the Chandra observations, we ex- 
tract high-resolution spectra corresponding to the dip and 
flaring intervals. We show the resulting spectra in Figure 
[2] As discussed in Section [2l we model the Chandra X- 
ray continuum with polynomials to account for calibration 
uncertainties. This technique effectively isolates narrow ab- 
sorption lines by characterizing the local continuum. We use 
Gaussian lines to measure the properties of the ionized iron 
absorption lines during each time interval; our results are 
shown in Table [2] 

Both Figure [2] and Table [2] reveal clear differences be- 
tween the accretion disc wind during the flare and dip inter- 



vals. From Figure[2] it can be seen that during the flares (top 
panel), we detect strong absorption lines from both Fexxvi 
(Lya: Is - 2p, 6.96 keV) and Fexxv (Hecv: Is 2 - ls2p, 6.7 
keV). The significance of the Fexxv and Fexxvi lines (~ 6a 
and ~ llcr) indicates a relatively strong, highly-ionized out- 
flow. During the dip (middle panel), Fexxvi is much weaker 
and Fe xxv is not significant at all, although it can be seen as 
a faint blip in the residuals (bottom panel). Although the 5- 
8 keV continuum flux decreases (Tabled from the flares to 
the dip, the bolometric luminosity remains roughly constant 
(Section HP)) . 

Interestingly, Table shows that aside from the reduc- 
tion in line flux, there are few if any quantitative differ- 
ences between the iron absorbers during these two parts of 
the fi cycle. During the flares, the wind is blueshifted by 
900+ 23 o km s _1 and has a line-of-sight velocity dispersion 
of cr v = 660+ 24 o km s _1 . In the spectrally-hard dip, the 
blueshift is 1340t^ km s 1 and the lines are unresolved, 
but we constrain the velocity dispersion to be < 1310 km 
s - . These results are consistent with the long -term behav- 
ior of the wind that we reported in iPaper J . Perhaps the 
most surprising result, however, is that despite the low sig- 
nificance of Fe xxv in the dip, the dip and flare spectra are 
statistically consistent with a constant Fexxv/Fexxvi ra- 
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Figure 2. The Chandra HETGS spectra and residuals for the 
flaring (top panel) and dip (middle panel) intervals of the state 
of GRS 1915+105. We detect strong absorption lines from Fexxvi 
during both time periods, and Fexxv during the X-ray flaring. 
These spectra include portions of observations 6579 and 6580, 
coadded to improve our sensitivity. The residuals are calculated 
with the line normalizations set to zero. 



tio. This suggests that the ionization state of the plasma 
in the wind might not vary from the dip to the flares, de- 
spite known changes in the hardness of the X-ray spectrum 
(Section E/JJ. 

We can substantiate this claim using a simple multi- 
plicative Gaussian model based on chapter 9 of lDrain el (|201ll ) 
that replaces the line flux parameter with the ion column 
density. First, we fit both phases of the cycle separately, al- 
lowing the Fexxvi/Fexxv ratio to vary between the dip and 
the flares. Then we fit them together, requiring both inter- 
vals to have the same Fexxvi/Fexxv ratio and line width. 
With the line ratios free, we achieve x 2 /f = 884.0/791 = 
1.12; with the line ratios and widths tied together, we obtain 
the same reduced x 2 = 885.8/793 = 1.12. In other words, 
there is no evidence (at the level of our line detections) that 
the ionization state of the wind changes between the dip and 
the flare. 

Using the results from our "tied" model, we find that 
the average Fexxvi/Fexxv ratio is r = 0.17too5; the 
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Figure 3. Chandra HETGS flux-resolved spectra and residuals 
for the flaring interval of the /3 state. As discussed in the text, 
the weakening of the Fexxv line can be taken as an indication 
that the ionization of the wind may increase slightly with the 
continuum flux during the flares. 

Fexxvi column density is A^F e 26 = 2.6tg;g x 10 lr cm -2 dur- 



ing the dip and A?Fe26 



8.0 



+ 1.6 
1.3 



x 10 cm during the 



flares. To convert from the ion column density to equivalent 
hydrogen columns, we as s ume elemental abundances from 
IWilms. Allen, fe McCravl |2000l ). Assuming only the ionic 
abundances of Fexxv and Fexxvi are non-zero, then these 
results imply an equivalent hydrogen column density in the 
wind of N H = TO ± 0.3 x 10 22 cm -2 during the dip and 
A H = 3.0l° g x 10 22 cm -2 during the flares. We obtained 
similar results with pho toionization models of the disc wind 
in the p state jPaper it . so we believe our simple measure of 
Ah here is plausible. We conclude that the column density 
of the wind decreases significantly from the flares to the dip, 
while its bulk dynamical and ionization properties remain 
constant. 

We can also ask whether or not the absorption lines vary 
within the flaring phase (our S/N is too low to perform a sim- 
ilar analysis of line variability in the dip). Using the Chandra 
lightcurve, we create new GTI files to divide the flare phase 
into three flux bands with roughly equal exposure times. The 
resulting spectra (see Figure [3| reveal significant differences 
in the wind between the lowest flux levels and the highest 
flux levels. Visible differences include apparent decreases in 
the equivalent widths and line widths, and a decrease in the 
strength of the Fe xxv line relative to Fe xxvi with increas- 
ing continuum flux. There is no evidence for variations in 
the velocity of the wind, but the changes in the line width 
and the Fexxv/Fexxvi ratio are marginally significant at 
90% confidence (Table [2]). Thus it appears that the ioniza- 
tion state or density of the wind may be somewhat sensitive 
to the luminosity on short time-scales. 

3.2 RXTE PCA 

With sensitive coverage out to 45 keV, RXTE PCA spec- 
tra can also provide valuable insight into variable ionization 
processes, so we extract PCA spectra for the dip and flar- 
ing phases with the same GTIs used for the HETGS (see 
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Figure 4. RXTE PCA spectra of the /3-state dip and flare, for 
comparison to Figure [2] The dip spectrum is notably harder than 
the spectrum of the flaring phase, but the 1 ev - 1 MeV lumi- 
nosities are comparable at 10 39 ergs s —1 . This suggests that the 
changes in the X-ray absorber may require variations in the den- 
sity or geometry of the wind (see Section [57TJ . 

Figure [TJ. Our purpose here is not to explore the variability 
of the accretion disc and corona, which has been done for 
the B state (e.g. iMigliari fe Bellonil [iool ; iRothstein et alj 
2005), but to capture the average shape of the ionizing 
spectrum. With typical wind densities around 10 12 cm~ 3 
(|Lee et alj|2003 : lUeda et alj|2009l ; fNeilsen et al.ll201ll ). the 
reco mbination time-scale is of the order of a few tenths of 
ms IIKallman et al . 2009). Thus our treatment of the aver- 
age spectrum is justified because the wind can effectively 
respond instantly to changes in the ionizing spectrum. 

We find that the average dip spectrum can be 
reproduced wi th a model consi sting of cold absorp- 
tion (tbabs; Wil ms et al.l [2000l) . a disc component 



(ezdiskbb; IZimmerman et all 2005h. a Comptonized com 



ponent (nthcomp; Zdziarsk i. Johnson, fe Magdziara 1 19961 : 



IZvcki. Done, fc Smitblll999l ), and an emission line in the 5-7 
keV range (egauss). The disc is cool (~ 0.75 keV) with a 
normalization (~ 800) that implies an inner radius around 
120 km. The photon index is Y ~ 2.5, and the electron tem- 
perature is ~ 12 keV; we tie the seed photon temperature to 
the disc temperature. During the flaring interval, we replace 
nthcomp with a high-energy cutoff (highecut) and simpl, a 
convolution model that takes a seed spe ctrum and scatters a 
fraction of the photons into a power law (Ste iner et al.l l2009). 
Our fits to the average flare spectrum imply a slightly hot- 
ter, smaller disc (T ~ 1.1 keV, R ~ 100 km), a steeper 
power law (T ~ 3.8) comprising ~ 60% of the photon flux, 
and an e-folding energy > 75 keV (we set the cutoff energy 
to zero). The fit parameters themselves must be interpreted 
cautiously given the strong variability in the flare and dip, 
but our models provide a good description of the spectra 
(X 2 /^ = 222.9/264; see Figure g]). 

Again, our primary interest is in the influence of the 
broadband continuum on the ionization of the wind, par- 
ticularly the strength of the Fe xxv absorption line. To that 



end, we integr ate our models over the Fe xxv photoionization 
cross-section (|Verner et al.|[l996l ) to determine the expected 
ionizing flux. In the flares, we find ~ 1.3 x 10 -3 ionizing pho- 
tons s _1 cm -2 ; during the dip, we find ~ 20% more ionizing 
photons (~ 1.5 x 10~ 3 photons s" 1 cm -2 ). Both spectra have 
0.001-1000 keV luminosities near 10 39 ergs s -1 . Since the lu- 
minosity and ionizing flux changes are modest, it is conceiv- 
able that the ionization state of the wind remains constant 
between the dip and the flares. However, the impact of the 
ionizing flux on the lines also depends on the changing den- 
sity and geometry of the wind, as well a s its chara cteristic 
distance from the continuum source (see IPaper ij and ref- 
erences therein). We explore these factors in more detail in 
Section 15.11 



4 THE 7 STATE 

In this section, we turn our attention to the variability of 
the accretion disc wind in the 7 state (ObsID 6581; Table 
[IJ. In Paper I, we reported the average RXTE continuum 
to be very soft, with ~ 82% of the 3-18 keV X-ray lumi- 
nosity emitted below 8.6 keV (Lx ~ 13 x 10 38 ergs s" 1 ). 
In the time-averaged Chandra HETGS spectrum, we de- 
tected an Fexxvi absorption line with an equivalent width 
of -21.9±l% eV and a blueshift of lOOOl^o km s" 1 . 

Although no signifi cant jet activity has be en reported 
in the 7, 8, or fi states (|Klein-Wolt et al.ll2002h . the strong 
red-noise variability (bottom panel of Figure[TJ may still pro- 
vide clues about the links between radiation and outflows in 
accreting X-ray binaries. In what follows, we discuss the vari- 
ability of the wind in this state, which is largely unexplored. 

The erratic variability here poses so me difficulties 
for g rating studies, since phase-folding (e.g. iNeilsen et alj 
1201 ll) and filtering on hardness-hardness-intensity diagrams 
( Soleri et al.l 120081 ) are not possible. These difficulties are 
compounded by the telemetry saturation discussed in Sec- 
tion [5] The simplest approach, dividing the observation into 
several time intervals, reveals no statistically-significant vari- 
ability in the wind, so we opt to extract flux-resolved spectra 
based on the count rate in the Chandra lightcurve. Because 
the exposure time and S/N are already relatively small, we 
use a sliding box (e.g. ISchulz fc Brandtll2002l ; INeilsen et alj 
1201 ll ) to create five spectra, each covering roughly one- 
third of the count rate distribution (with overlap). While 
this method alleviates the difficulty in selecting physically- 
motivated GTIs, it does create some overlap between consec- 
utive spectra. We only use statistically-independent spectra 
for our variability calculations. 

We show the three statistically-independent spectra 
(71, 73, and 75) in Figure [S] all three spectra contain a 
strong Fexxvi line with a blueshift of 1000 — 1500 km s -1 . 
In contrast to the /3-state spectra, here there is little evi- 
dence for absorption by Fe xxv (we ignore the weak features 
near 6.6-6.7 keV in 73 and 75 because they do not show up in 
both of the HEG ±1 orders). The absence of Fexxv is plau- 
sible from an ionization perspective, given the high contin- 
uum flux in the 7 state. Our multiplicative Gaussian model 
( Section |3T|) indicates that Nf»26 ~6x 10 17 cm ~ 2 . If all 
the iron is locked up in Fe xxvi (again, plausible because we 
do not observe any other lines), this value of 2VF e 26 implies 
an equivalent hydrogen column density of Ah ~ 2 x 10 22 



cm . However, this hydrogen column should be regarded 
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Figure 5. The Chandra HETGS flux-resolved spectra and resid- 
uals for three distinct intervals in the 7 state. The continuum 
level changes, but the lines are very similar during these intervals. 
However, the complicated relationship between the absorbed line 
flux and the continuum (Section[4j indicates that photoionization 
effects are still important at these high luminosities. 

as a lower limit, since there could be a significant popula- 
tion of fully-ionized iron (which would increase the total iron 
column density without affecting the lines). 

As for the variability of the wind in the 7 state, it is 
clear from Figure [5] that the continuum flux is variable, but 
the lines themselves are difficult to distinguish because they 
have similar significance, absorbed fluxes, and blueshifts. 
From Table [2] however, the Fe xxvi line flux is highest at 
intermediate luminosities, and the equivalent width drops at 
the highest luminosities. We rule out a constant line flux at 
the ~ 97% confidence level. The variation of the equivalent 
width between all three spectra is not particularly signifi- 
cant, but the decrease in Wo in 75 relative to 71 and 73 has 
a formal significance of ~ 97%. In general, the lines do not 
appear to track the continuum flux directly, and so may be 
subject to other influences (Section [5J. 

5 DISCUSSION 

In GRS f9f5+105, the accretion disc wind is now known to 
vary on time-scales ranging from a few seconds to (nearly) 
decades. This variability appears particularly in the column 
density and ionization state of the wind, but is also evident 
in the gas density, velocity, and velocity dispersion. In the 
rest of this section, we discuss the implications of our results 
for the complex relationship between the properties of the 
wind and the luminosity of the central X-ray source, which 
is responsibl e not only for driving the wind, but also for ion- 
izing it (e.g. iBeeelman et al.lll983l : IProea fe Kallmanll2002l : 
iNeilsen et alj|201ll ). 



5.1 The (3 state: dip vs. flare 

First we consider the dip and flare intervals of the /3 state 
described in Section T3. II To recap briefly, we find that when 
moving from the flaring phase into the hard X-ray dip, the 
column density of the wind decreases significantly (~ 3x), 



while its velocity and velocity dispersion remain constant. 
Furthermore, the ionization state of the absorber does not 
appear to change despite changes in the shape of the ionizing 
spectrum. We argued in Section ^. 2l that changes in the den- 
sity or geometry of the wind might therefore be important 
in determining the properties of the absorber. 

For a toy model of this variability, suppose that the 
wind is launched from the disc as a sphere expanding with 
constant velocity v, and suppose further that the wind is 
launched solely during the bright, soft, X-ray flares. In other 
words, no wind is driven off the disc during the hard X-ray 
dip. For simplicity, we assume the wind is spherical, but as 
long as its solid angle is roughly constant this toy model will 
still apply. By the end of the flaring phase of the f3 cycle, 
which we define as T = 0, this sphere has a characteristic 
radius Ro that is comparable to the typical launch radius of 
the wind (see the top panel of Figure [5]). During the course 
of a flaring interval, the accretion disc launches N particles 
into the wind, so that the expected column density is 

During the dip, an average time AT later, when the charac- 
teristic radius of the wind is 



R!=Ro + vAT, 
the expected column density is 
N 



AtvR/ 



N 



(2) 

(3) 
(4) 



4tt(R +vAT) 2 ' 

Given the ratio of column densities in and out of the hard 
X-ray dip, as well as the wind speed and AT, we can solve 
for the wind launch radius Ro : 

'AT 



Ro 



4> 



(5) 



where "I> = y Nno /Nm. Using the values reported in Section 
13.11 Table [2] and an average mid-flare to mid-dip AT = 626 
s, we calculate Ro = (8 ± 3) x 10 10 cm. The error comes 
from propagating la uncertainties for each of the measured 
quantities through Equation ((5JI . The implied gas density in 



, and the expected 
\ This is 



the wind is n e = (1.4 ± 0.5) x 10 cm" 
ionization parameter is around 10 5 ' 2 ergs cm s 
rather high for the observed lines, but could be reduced if 
the wind is clumpy or accelerates in the dip, as hinted at in 
Table The average mass loss rate in the typical ~ 800 s 
flare phase is M w i n d = (7 ±5) x 10 18 g s _1 ; the instantaneous 
mass loss rate from the disc could be significantly higher if 
the wind is only launched during a fraction of the flare phase, 
or if the ionic abundance of completely-ionized iron is not 
negligible during the dip. 

Of course, the fact that this toy model returns sensible 
results does not imply or require that it is correct. We could 
alternately suppose, for example, that the wind is produced 
at a constant rate throughout the cycle. In this case, we can 
replace $ in equation 5 with \ = & * yl + ATdi p /ATfl are . 
Here the ATs represent the duration of the dip and flare 
intervals. For our data, $ ~ 1.7 and x ~ 2.6, so this mod- 
ification to our model reduces Ro by a factor of ~ 2. The 
density n e increases by a similar factor, but because the ion- 
ization parameter £ scales with R~ 2 , the modified model 
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Figure 6. Toy model for the evolution of the column density 
in the state of GRS 1915+105. (TOP): the flaring interval, 
where the (spherical) wind has characteristic radius Ro and a 
high column density (indicated by dark shading). (BOTTOM): 
during the dip, a time AT later, the wind has expanded and 
its column density has decreased accordingly (lighter shading). 
In Section 15.11 we use this model and our measurements of the 
column density to estimate Ro ~ 8 X 10 10 cm. 



implies even more ionization than our original model (where 
the wind shuts off in the dip). Overall, both are reason- 
able, but the excess ionization in the modified (constant 
wind) model is potentially problematic. Again, the ioniza- 
tion parameter could be reduced if the wind is clumpy, 
but the modified model requires an even clumpier wind 
than our original model, since the ionization parameters are 
^modified ~ 2£ or iginai- Furthermore, we know from previous 
work that hard X-ray spectra (similar to that observed in the 
dip) tend to indicate activity in jets rather than disc winds 
(e.g. iKlein-Wolt et all 12001 iFender. Belloni. fc Gallol 12004 
iNeilsen fe Leeh o09 and references therein). This would lead 
us to the prediction that the hard dip would be more likely 
to produce a jet than a wind. For these reasons, we suggest 
that our original model for the wind variability in the j3 state 
is preferable, and that the wind actually does turn off during 
the dip. 

As in our previous work on the persistent properties 
of the wind on long time-scales and its variability on short 
time-scales dNeilsen fc Le e 2009 ; Neilsen et al.ll201ll ; see also 
iLuketic et alJl201C| ). our results are consistent with a wind 
driven primarily via a combination of Compton heating of 
the outer accretion disc and radiation pressure on free elec- 
trons. Furthermore, our results hint at a unifying process 



driving the wind. In the heartbeat (p) state, which consists 
of bright X-ray pulses every ~ 50 seconds, we found that 
each X-ray pulse heats up a large re gion of the disc an d 
drives new gas into the line of sig ht |Neilsen et ai"1l201ll ). 
Since the flares in the /3 state are similar to the pulses in the 
heartbeat state (as will be demonstrated definitively in our 
forthcoming analysis of a new Chandra, RXTE, EVLA and 
Gemini campaign; Neilsen et al. 2012, in preparation), the 
mechanisms producing the winds are likely similar as well. 

Yet the wind in the ft state also exhibits a critical dif- 
ference from the wind in the heartbeat state, for our toy 
model of the /3 state suggests that the wind actually shuts 
off during the hard X-ray dip. We find no evidence in the 
line ratios for over-ionization of the absorber during the dip, 
and the ionizing flux is comparable during both phases of the 
cycle. It seems improbable that such modest changes in the 
ionizing and bolometric luminosities could lead to enough 
over-ionization to reduce the apparent wind column density 
so significantly. Thus we conclude that there is a long in- 
terval where no wind is launched from the disc, a fact that 
is rather counterintuitive because the wind absorption lines 
never fade completely. But during the dip, we actually de- 
tect the remnant of the wind launched by prior flares (in the 
form of significantly weaker absorption lines). 

This result poses an intriguing puzzle: why doesn't the 
disc produce any wind during the X-ray dip? Our data do 
not provide any clear answers to this question, but the dip 
and the flares have similar average X-ray luminosities, so it 
seems improbable that the wind is quenched solely due to 
changes in the heating or energetics of the outer disc. 

One possibility is that the highly-ionized wind in GRS 
1915+105 is primarily produced by short but intense pulses 
of radiation, as observed in the heartbeat and /3 states. If this 
impulsive driving process launches partially- or temporarily- 
ballistic parcels of gas into the line of sight, it might explain 
how the Fexxvi/Fexxv ratio can rem ain roughly constant 
over a large range in luminosity (see iKallman fc Bautistal 
2001): ionization would only become important as the wind 
rises well abo ve the disc. Along th is line of reasoning, it is in- 
teresting that lMiller et al. (2006b) detected a wind with very 
similar ionization properties in H1743-322 while the source 
was undergoing moderate ~ 300 s oscillations. Still, a general 
link between wind formation, ionization, and X-ray variabil- 
ity leaves much to be explained. For example, if variability 
is important, why are these winds rarely observed during 
X-ray hard states, where the rms X-ray variability is high? 
What determines whether the wind is highly ionized, as in 
H1743-322 and the /3, 7, and p states of GRS 1915+105, o r 
cooler with ma ny lines, as in GX 13+1 (|Ueda et al.ll2004T ). 
GRO J165 5-40 jMiller et al.ll2006al ). and the <j> state of GRS 
1915+105 (|Ueda et al.ll2009l )? 

With these questions unanswered, a simpler explanation 
for the absence of a wind is that the outer accretion disc is 
shadowed during the hard X-ray dip and does not receive 
the full irradiating force of the X-ray luminosity. This could 
be possible if, for example, the inner disc is puffed up dur- 
ing the dip, so that low inclinations are shielded from the 
X-ray emission (see also lPonti et al.ll201lf l. In this explana- 
tion we find an exciting link to the formation of jets, since it 
is believed that steady jets may require geometrically-thick 
discs (e.g. M eieil [200 ll ). To be clear, t he /3 state produces 
discrete synchrotron blobs (baby jets; I Fender et al.l 1 19971 : 
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iPoolev fc Fended Il997l ; iMirabel et all 1 19981 ), not the com- 
pact, flat-spectrum jets for which the thick disk theory was 
developed. However, a thick disk may also be required to pro- 
duce other kinds of jets, like those observed after the hard dip 
in the /3 state. Whether th e baby jet is launched co ntinuously 
during the hard dip (e.g. iKlein-Wolt et al. 20021') or at th e 



moment of the "trigger spike" (e.g. Rothstein et al.l f2005). 
it seems reasonable to assume that the dip in the lightcurve 
coincides with changes in the structure of the inner disc that 
ultimately lead to the ejection of a synchrotron blob. These 
same changes may therefore hide a significant fraction of 
the central X-ray luminosity from the wind-producing outer 
disc. 

As for the relative timing of the disc and wind 
variability, we note that the ft state is associated with 
periodic changes in t he disc mass accretion rate (e.g. 
iMigliari fc Bellonil l2003 ) , whether or not their origin is com- 
pletely understood. It is expected that mass accretion rate 
and scale height in the inner disc can evolve on viscous 
and t hermal time-scales, respectively |Frank. King, fe Raind 
2002). The thermal time-scale is comparable to the dynam- 
ical time-scale, so the scale height of the disc can change 
rapidly relative to the accretion rate at any radius. To 
be quantitative, we consider the radius i? v i sc at which the 
viscous time t v i sc is equal to the j3 state's 30-minute cy- 
cle (i.e. we assume that the j3 state is in part a viscous 
instability like many other states of GRS 1915+105; e.g. 
Navakshin. Rappaport. fc Meliall2000l ). Using Equation 1 of 
Belloni et al.l l| 19971 ). we find that R viac < 7 x 10 9 cm even if 



M is as high as 10 g s (assuming a viscosity parameter 
a = 0.01 and black hole mass M = 14 M©). At this ra- 
dius, the thermal time-scale is ~ 0.1 second. Thus any part 
of the disc participating in 30-minute viscous variations can 
change its scale height almost instantly, and in any case on 
time-scales much shorter than the duration of the X-ray dip. 
Furthermore, we know that winds driven by thermal and ra- 
diation pressure can res pond to c hanges in the luminosity 
in as few as five seconds l|Paper lit) . Therefore, we conclude 
that our shadowing explanation is plausible. 

If this scenario is correct, it would explain the uncanny 
ability of the wind to "know" about the state of the inner 
disc, i.e. why winds might be preferentially formed when jets 
are absent. Because the radiation from the disc launches the 
wind, the wind can effectively respond to changes in the 
driving luminosity as they occur. There may be delays due 
to light travel time, but this is only a few seconds and is 
short compared to the duration of the dip. In short, the 
remarkable quenching of the wind during the hard X-ray 
dip is a clear consequence of the radiative links between the 
inner and outer parts of the accretion flow, consistent with 
our results for the he artbeat state on time-scales of 5 seconds 
(|Neilsen et al.ll201ll ). 

5.2 Flux-resolved studies 

The complex radiative links between the wind and the inner 
and outer accretion disc are also apparent in the flux depen- 
dence of the properties of the wind absorption lines. For the 
P and 7 states, we summarize our results in Figure [T) Note 
that any two points in Figure that have the same symbol 
and fill style are statistically independent. At this level, we 
do not have the signal to statistically distinguish changes in 
the wind speed, but this figure clearly demonstrates that the 
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Figure 7. Flux-resolved properties of the accretion disc wind in 
the and 7 states of GRS 1915+105. The /3 state points are 
marked with circles, while the results for the 7 state are marked 
with squares. For a given state, symbols with the same fill style 
are statistically independent. Top panel: the Fexxvi line flux, 
which generally increases with increasing continuum flux. Middle 
panel: the wind blucshift, which is consistent with a constant over 
all flux. Bottom panel: the Fexxvi equivalent width, which is 
mostly constant; see Section \E\ for a discussion of deviations from 
a constant equivalent width. 



line flux and equivalent width exhibit a dependence on the 
observed continuum flux that is non-linear and not necessar- 
ily even monotonic. 

Part of this non-linearity is due to the behavior dis- 
cussed in Section l5,ll in which residual absorption lines from 
a wind launched in some previous epoch may still be present 
in the X-ray spectrum. As discussed in the preceding section, 
wind driving may be sensitive to both the X-ray luminosity 
and the state and structure of the inner accretion flow. 

For the most part, the line flux rises smoothly and the 
equivalent width is remarkably constant over a factor of 
nearly five in continuum flux. This would seem to indicate 
that the Fexxvi column density changes very slowly with 
luminosity, even across different states. A direct comparison 
between (3 and 7 is difficult because of possible changes in 
line saturation, and because the lines are generally resolved 
in the 7 state but generally not in the /3 state. Nevertheless, 
the smooth evolution of the wind properties suggests a com- 
mon origin (i.e. radiative/thermal driving), in which case we 
can expect the accretion disc wind and the continuum flux 
to be strongly linked. 
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For example, the weakness or absence of Fe xxv in the 7 
spectra (relative to (3) does suggest that as the X-ray lumi- 
nosity increases, the wind does become progressively more 
ionized. This conclusion is supported by the behavior of 
the Fexxv/Fexxvi ratio in our flux-resolved ft flare spec- 
tra. In addition, because the abundance of Fexxvi should 
drop sharply due to over-ionization above some critical lu- 
minosity, we believe the drop in the Fexxvi line flux and 
equivalent width after 73 (Figure [7]) further reinforce the 
conclusion that the X-ray luminosity is a major determining 
factor of the wind properties. 

On the other hand, it remains clear that there is more to 
disc winds than the luminosity that drives them. For exam- 
ple, the density of the wind in GRS 1915+105 se ems to be 
independent of luminosity at around 10 12 c m" 3 (|Lee et alj 
12003 ; lUeda et al l 120091 ; iNeilsen et all 1 20 lit this work). Or 
consider that the equivalent width of the iron line is highest 
in Pi, which has the lowest continuum flux studied in this 
work. It may be that different driving mechanis ms are impor- 
tant at different luminosi ties, i.e. MHD winds (|Progall2000l ; 
iMiller et ail l2006al . l2008f ). or it may be that flux-resolved 
studies like the one presented here fail to capture the true 
time-dependence of variable winds. In the heartbeat state, 
the wind is launched by high luminosities but doesn't rise 
into o ur line of sight until after the X-ray luminosity has 
fallen (|Paper III) ; the same may also be true during the /? 
flares and the 7 state, whose variability is even more erratic. 



6 CONCLUSIONS 

In this paper, we have presented a time- and flux-resolved 
high-resolution spectral analysis of GRS 1915+105 in the /3 
and 7 states. Although only the j3 state is associated with 
transient jet formation, both states exhibit strong erratic 
X-ray variability, which we have used to probe the relation- 
ship between the central X-ray luminosity and the proper- 
ties of the accretion disc wind. In particular, our new results 
confirm our previous finding that radiation mediates rapid 
interactions between the inner accretion disc and the outer 
accretion disc. This interaction is most evident in the f3 state, 
where the wind is actually quenched during the hard X-ray 
dip, possibly because the irradiation of the outer disc is af- 
fected by changes in the structure of the inner disc leading 
to the formation of a transient jet. To the best of our knowl- 
edge, this is the first time that a wind has been observed to 
turn off within a single observation, and the mechanism (M 
and/or structural changes in the inner disc) requires intri- 
cate links between radiation, winds, and jets. 

This result is therefore an exciting step towards a com- 
plete understanding of wind-jet interactions on all time- 
scales. With our sensitivity and limited time resolution, our 
best estimates place the wind in the /3 state ~ 10 11 cm from 
the black hole, leaving it effectively unable to affect the tran- 
sient jet on such short time-scales. However, if future mul- 
tiwavelength studies unaffected by telemetry saturation are 
able to place better constraints on the location of the wind, 
and require it to be very close to the black hole, then the 
P state might produce the first evidence for direct wind-jet 
interactions on any time scale. Such a result would signifi- 
cantly enhance our understanding of accretion and ejection 
as profoundly inter-related processes around black holes. 
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